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Abstract A search for neutrinoless (3/3 decay processes 
accompanied with Majoron emission has been performed 
using data collected during Phase I of the GERmanium 
Detector Array (Gerda) experiment at the Laboratori 
Nazionali del Gran Sasso of INFN (Italy). Processes 
with spectral indices n = 1, 2, 3, 7 were searched for. 
No signals were found and lower limits of the order of 
10 23 yr on their half-lives were derived, yielding sub¬ 
stantially improved results compared to previous ex¬ 
periments with 76 Ge. A new result for the half-life of 
the neutrino-accompanied (3f3 decay of 76 Ge with sig¬ 
nificantly reduced uncertainties is also given, resulting 
in Ty 2 = (1.926 ± 0.095) • 10 21 yr. 

Keywords double beta decay • Majoron emission • 
enriched 76 Ge 

PACS 23.40.-s /3 decay; double (3 decay; electron 
and muon capture • 14.80.Va majorons • 21.10.Tg 
Lifetimes, widths • 27.50.+e mass 59 < A < 89 

1 Introduction 

Neutrinoless double beta (Cfo/3/3) decay is regarded as 
the gold-plated process for probing the fundamental 
character of neutrinos. Observation of this process would 
imply total lepton number violation by two units and 
that neutrinos have a Majorana mass component. Al¬ 
though the main focus of experimental efforts lies on the 
detection of Ovf3(3 decay mediated by light Majorana 
neutrino exchange, there are also many other proposed 
mechanisms which are being searched for. Some exotic 
models predict Cfo/?/3 decays proceeding through the 
emission of a massless Goldstone boson, called Majoron. 
Predictions of different models depend on its transfor¬ 
mation properties under weak isospin, singlet [1], dou¬ 
blet [2] and triplet [3] . Precise measurements of the in¬ 
visible width of the Z boson at LEP [4] greatly disfavour 
triplet and pure doublet models. Several new Majoron 
models have been developed subsequently in which the 
Majoron carries leptonic charge and cannot be a Gold- 
stone boson [5,6] or in which the 0vf3(3 decay proceeds 
through the emission of two Majorons [7]. 

All these models predict different shapes of the two 
emitted electrons’ summed energy spectrum. The pre¬ 
dicted spectral shapes are essentially defined by the 
phase space of the emitted particles: 

^-G-iQpp-K)" ( 1 ) 
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where K is the summed energy of the two electrons, G 
is the phase space, Q/3/3 is the Q value of the 0vf3[3 decay 
and n is the spectral index of the model. Single Majoron 
emitting (3f3 decays can be roughly divided into three 
classes, n = 1, n = 2, and n = 3. Double Majoron 
emitting decays can have either n = 3 or n = 7. Their 
characteristic spectral shapes differ from that of two- 
neutrino /3/3 decay (2uf3f3), for which n = 5. This allows 
for discrimination between the processes. 

Experimental searches for /3/3 decay mediated by 
emission of one or two Majorons (0fo?/?x) have been 
performed by the Heidelberg-Moscow experiment (HdM) 
for 76 Ge [8,9]; by Nemo-2 and Nemo-3 for 100 Mo, 116 Cd, 
82 Se, 96 Zr, 130 Te [10,11]; by ELEGANT V for 100 Mo [12]; 
by DAMA [13] and by KAMLAND-Zen for 136 Xe [14], 
None of these experiments have seen an excess of events 
that could be interpreted as a Majoron signal; they re¬ 
ported lower limits on the half-lives of the processes 
that involve Majoron emission. 

The 2vf3f3 decay process conserves lepton number 
and is independent of the nature of the neutrino. It has 
been detected for eleven nuclides so far, with measured 
half-lives (TfY 2 ) in the range of 7 x 10 18 — 2 x 10 24 yr [15, 
16,17]. The knowledge of Ty 2 allows for extraction of 
the nuclear matrix element, A4 2y , which can provide 
some constraints on that of Ozx/3/5 decay, A4° 1 ', if the 
evaluations of A4 for the two processes are performed 
within the same model [18,19]. 

This paper reports on the search for neutrinoless 
double beta decay of 76 Ge with Majoron emission ( 0vf3(3x ) 
and a new analysis of the half-life of the 2v(3(3 decay of 
76 Ge using data collected by the Gerda experiment 
during its Phase I. 2zx/3/5 decay is a well established and 
previously observed process, while Ovfifix decay is a hy¬ 
pothetical one. In the first case the half-life is extracted, 
while for the second one a limit is set. This leads to 
slightly different approaches in the analyses leading to 
different data sets and background components being 
used. 


2 The GERDA experiment 

The main aim of the Gerda experiment [20] at the 
Laboratori Nazionali del Gran Sasso (LNGS) of INFN 
in Italy is to search for Ozz/3/3 decay of ' 6 Ge. The core of 
the setup is an array of high-purity germanium (HPGe) 
detectors made from isotopically modified material with 
76 Ge enriched to ~86% ( e " r Ge), mounted in low-mass 
copper supports (holders) and immersed in a 64 m 3 
cryostat filled with liquid argon (LAr). The LAr serves 
as cooling medium and shield against external back¬ 
grounds. The shielding is complemented by water in a 
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tank of 10 m in diameter which is instrumented with 
photomultipliers to detect Cherenkov light generated in 
muon-induced showers [20]. 

The array of HPGe detectors is arranged in strings. 
Each string is enclosed with a cylinder, made from 60 pm 
thick Cu foil, called mini-shroud, to mitigate the back¬ 
ground coming from the decay of 42 Ar present in the 
LAr. Moreover, in order to prevent contamination from 
radon within the cryostat, a cylinder, made from 30 pm 
thick Cu foil, called radon-shroud, separates the central 
part of the cryostat, where the detectors are located, 
from the rest. The HPGe detector signals are read out 
with custom-made charge sensitive preamplifiers opti¬ 
mized for low radioactivity, which are operated close to 
the detectors in the LAr. The analog signals are digi¬ 
tized with 100 MHz Flash ADCs (FADC) and analyzed 
offline. If one of the detectors has an energy deposition 
above the trigger threshold (40-100 keV), all channels 
are read out. Reprocessed p-type coaxial detectors from 
the HdM [21] and Igex [22] experiments were operated 
together with Broad Energy Germanium (BEGe) type 
detectors manufactured by Canberra [23,24]. 

As explained in section 5, some background compo¬ 
nents have different effects on the two detector types 
due to their peculiar geometry. A schematic drawing 
of a coaxial detector type is shown in the top part of 
Fig. 1, while the lower part depicts that for a BEGe 
type detector. 


60-80 mm 



Fig. 1 Schematic sketch of a coaxial HPGe detector (top) 
and a BEGe detector (bottom) with their different sur¬ 
faces and dead layers (drawings not to scale), adapted from 
Ref. [25]. 


3 Data taking and data selection 

Phase I data taking lasted from November 9, 2011, to 
May 21, 2013. The total exposure collected comprises 
19.2 kg-yr for the coaxial detectors and 2.4 kg-yr for 
the BEGe detectors. In this paper, the entire expo¬ 
sure collected by the BEGe detectors (BEGe data set) 
and 17.9 kg-yr from the coaxial detectors (golden data 
set) are used [25,26]. For the coaxial detectors, a data 
set collected for 1.3 kg-yr exposure during a restricted 
time period around the deployment of the BEGe detec¬ 
tors is discarded due to a higher background level. Also 
one of the coaxial detectors, RG2, is not considered for 
the data analysis starting from March 2013, as its high 
voltage had to be reduced below depletion voltage due 
to increased leakage current. The energy calibration of 
the detectors was performed using the information from 
dedicated calibration runs. For these calibration runs, 
three 228 Th sources were lowered to the vicinity of the 
detectors. The stability of the energy scale was moni¬ 
tored by performing such calibration runs every one or 
two weeks. Moreover, the stability of the system was 
continuously monitored by injecting charge pulses into 
the test input of the preamplifiers. Using physics data, 
the interpolated FWHM values at Qgg averaged with 
the exposure are (4.8 ± 0.2) keV for the coaxial detec¬ 
tors and (3.2 ± 0.2) keV for the BEGe detectors. 

All steps of the offline processing of the Gerda data 
were performed within the software framework Gela- 
tio [27]. The energy deposited in each detector was 
extracted from the respective charge pulse by applying 
a approximate Gaussian filter [28]. Non-physical events, 
such as discharges, cross-talk and pick-up noise events, 
were rejected by quality cuts based on the time position 
of the rising edge, the information from the Gaussian 
filter, the rise time and the charge pulse height, which 
must not exceed the dynamic range of the FADCs. Pile- 
up and accidental coincidences were removed from the 
data set using cuts based on the baseline slope, the 
number of triggers and the position of the rising edge. 
The rate of pile-up and accidental coincidence events 
is negligible in the Gerda data due to the extremely 
low event rate. The loss due to mis-classification by the 
quality cuts was <0.1 % for events with energies above 
1 MeV. All events that come within 8 ps of a signal 
from the muon veto were rejected. Finally, only events 
that survive the detector anti-coincidence cut were con¬ 
sidered. This means, that all events with an energy de¬ 
position > 50 keV in more than one detector in the ar¬ 
ray were not taken into account. Since 2i//3/3 and 0u/3f3x 
events release their energy within a small volume inside 
the detectors, almost no signal events were lost by this 
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cut, while a part of the 7 -induced background events 
were rejected. 


4 Analysis Strategy 

The two analyses described in this paper are different 
in the sense that for 2 zz/3/3 decay a parameter is ex¬ 
tracted for a well established and known process, while 
in the case of the search for Qvftflx decay limits for 
a hypothetical process are set. In order to minimize 
the systematic uncertainties for the extraction of the 
T.jy 2 it is favorable to use a well defined and controlled 
subset of the data and to use only well identified back¬ 
ground processes. For 0^/3 fix limit setting it is favorable 
to maximize the exposure and to take into account all 
known possible background processes that can not be 
unambiguously detected but could mimic Qvfifix decay. 

For the Ty 2 analysis the golden data set (17.9 kg-yr) 
with the coaxial detectors is used in order to have a 
large data sample obtained in well controlled experi¬ 
mental conditions. The Majoron analysis uses both the 
golden data set and the BEGe data set for a total expo¬ 
sure of 20.3 kg-yr in order to maximize the sensitivity. 

The background model for the T ^ 2 analysis uses a 
minimal number of components, assuming all sources 
near to the detectors [25,29]. For the Majoron analysis, 
an expanded model is used [30] , taking into account also 
additional medium and far distant positions for some of 
the sources. This becomes necessary when searching for 
rare processes such as Majoron emission, where all pos¬ 
sible sources of background which could simulate the 
exotic process have to be considered. Therefore, even 
the slight differences resulting, for example from a vari¬ 
ation of the source position, have to be evaluated. 

In both analyses, the experimental spectra of the 
coaxial and BEGe detectors are analyzed using the Baye¬ 
sian Analysis Toolkit (Bat) [31]. 


5 The background model 

The background sources considered in the models were 
identified by their prominent structures in the energy 
spectra and were also expected on the basis of ma¬ 
terial screening measurements. The spectral shapes of 
individual background contributions were obtained by 
using a detailed implementation of the experimental 
setup in the Monte Carlo (MC) simulation framework 
MaGe [32], A Bayesian spectral fit of the measured 
energy spectrum with the simulated spectra was per¬ 
formed in an energy range from 570 keV up to the end 
of the dynamic range at 7500 keV. The low energy limit 


is motivated by the /3-decay of 39 Ar, which gives a large 
contribution up to its Q/ 3 -value of 565 keV. 

The following background components were used for 
the extraction of the Tfj. 2 (minimum model in Refs. [25, 
29]): (1) 76 Ge 2vfifi decay, (2) 214 Bi, 228 Ac, 228 Th, 60 Co 
and 40 K decays in the close vicinity of the detectors 
(<2 cm, represented by decays in the detector holders 
in the MC simulation), (3) decays of 60 Co inside the 
detectors, constrained by the maximum expected ac¬ 
tivity from their cosmogenic activation history, (4) 42 K 
decays in LAr assuming a uniform distribution, (5) a- 
model that accounts for a decays originating from 210 Po 
and 226 Ra contaminations on the p + surface of the de¬ 
tectors as well as from 222 Rn in the LAr, and finally 
( 6 ) 214 Bi decays on the p + surface, constrained by the 
estimated 226 Ra activity from the a-model. 

The parameters of all components besides the con¬ 
strained ones were given a flat prior probability dis¬ 
tribution. There are no strong correlations between the 
model parameters since all considered background com¬ 
ponents have characteristic features such as 7 -ray lines 
or peak-like structures at different energies. The ratios 
of the 7 -ray line intensities from the individual consid¬ 
ered background sources suggest contaminations dom¬ 
inantly in locations close to the detectors. Hence, the 
minimum model takes into account only the close-by 
source locations. Nevertheless, the screening measure¬ 
ments indicate contaminations of materials in farther 
locations as well. An additional contribution can come 
from 42 K decays at or near the detector n + surfaces 
(see Fig. 1) with a specific activity higher than that for 
the uniform distribution assumption. This component 
is the dominating one for the BEGe data set, as the 
thinner dead layer thickness of BEGes of roughly 1 mm 
allows penetration of the electrons emitted in the decay 
of 42 K to the active volume, while for coaxial detectors 
the dead layer thickness of ~2 mm efficiently shields 
this background component. 

The spectral shapes of the contributions from the 
background sources without significant multiple 7 peaks 
at different source locations differ only marginally. This 
makes it impossible to pinpoint the exact source loca¬ 
tions given the available statistics of the measured spec¬ 
tra. Therefore, variations of the source locations for the 
considered decays were taken into account when evalu¬ 
ating the systematic uncertainty on T 2 )^. 

For the Majoron analysis additional background com¬ 
ponents were used [30], including also medium and far 
distant contributions. For the coaxial detectors 42 K on 
the n + and on the p + contacts was added to the list of 
the close sources of the previous background model. For 
medium distances, i.e. between 2 cm and 50 cm from 
the detectors, contributions from the following sources 



5 


were added: 214 Bi, 228 Tlr and 228 Ac. A 228 Th contam¬ 
ination was chosen as a representative for far distant 
sources (above 50 cm). Whenever possible, screening 
measurements were used to constrain the lower limit of 
the expected background events. 

In the Majoron analysis, also the data collected with 
the BEGe diodes were used in order to maximize the 
exposure. Consequently, the background model devel¬ 
oped for these detectors was used [25,30]. The same 
close, medium and far distant sources as for the coaxial 
detectors were used. 68 Ge was added as internal source. 
This was necessary in order to take into account the 
cosmic activation of the germanium due to the recent 
production of these diodes. 


6 Determination of the half-life of 2v>f3f3 decay 


6.1 Analysis 


TheT 2 ; 2 of 2vfi ft decay of 76 Ge was determined consid¬ 
ering the golden data set of Phase I, amounting to an 
exposure of 17.9 kg-yr, and using the background model 
prediction for the contribution of the 2v/3/3 spectrum to 
the overall energy spectrum. Details of the background 
analysis can be found in Ref. [29]. 

The global fit for the background modeling was per¬ 
formed on the summed energy spectrum of the coaxial 
detectors using a bin width of 30 keV. Thus, the scal¬ 
ing parameter of the 2V/3/3 spectrum in the model, Wf*, 
gives the number of events in the 2u/3/3 spectrum in the 
fit window of 570-7500 keV for all detectors. Using this 
result for the number of measured 2vfifi events, the 
half-life is calculated as 


N de 


rpiv _ (In 2) Na \ i\j + f -fit 

^1/2 jvrfit M,t U 7 6,* L JAV,1 £ AV,i 

rn e nr tv 2 „ , =1 

+ (1 — /AV.i) ^DL,i] 


( 2 ) 


where Na is Avogadro’s constant and m en r = 75.6 g is 
the molar mass of the enriched material. The summa¬ 
tion runs over all the detectors ( Ndet ) considered in the 
data set. All detector related parameters like the de¬ 
tector mass the time of the data taking for each 

detector (tf), the fraction of 76 Ge atoms (/re,*), the ac¬ 
tive volume fraction (/av.z), an d the detection efficien¬ 
cies in the active volume (e®Vi) an d in the dead layer 
( £ dl i) are taken into account separately for the indi¬ 
vidual detectors. All values are listed in Table 1. The 
efficiency e^ Vi (sjf L fi) corresponds to the probability 
that a 2vj3fi decay taking place in the active volume 
(dead layer) of the detector deposits detectable energy 
in the fit window considered for the background model. 


Table 1 Parameters for the coaxial detectors (upper part) 
and for the BEGe detectors (lower part): live time, t, total 
mass, M, the fraction of 76 Ge atoms, / 76 , and the active 
volume fraction, f av■ For the coaxial detectors, the first un¬ 
certainty on fact is the uncorrelated part, the second one the 
correlated contribution. The values for M, fyg and /_av are 
taken from Ref. [25]. 


detectors 

t 

M 

h 6 

Sav 


[days] 

[kg] 

[%] 

[%] 

enriched coaxial detectors 

ANG2 

485.5 

2.833 

86.6 ±2.5 

87.1 ±4.3 ±2.8 

ANG3 

485.5 

2.391 

88.3 ± 2.6 

86.6 ±4.9 ±2.8 

ANG4 

485.5 

2.372 

86.3 ± 1.3 

90.1 ±4.9 ±2.9 

ANG5 

485.5 

2.746 

85.6 ± 1.3 

83.1 ±4.0 ±2.7 

RG1 

485.5 

2.110 

85.5 ± 1.5 

90.4 ±5.2 ±2.9 

RG2 

384.8 

2.166 

85.5 ± 1.5 

83.1 ±4.6 ±2.7 

enriched BEGe detectors 

GD32B 

280.0 

0.717 

87.7 ± 1.3 

89.0 ± 2.7 

GD32C 

304.6 

0.743 

87.7 ± 1.3 

91.1 ±3.0 

GD32D 

282.7 

0.723 

87.7 ± 1.3 

92.3 ± 2.6 

GD35B 

301.2 

0.812 

87.7 ± 1.3 

91.4 ±2.9 


The detection efficiencies, on average = 0.667 and 
= 0.011, are obtained through dedicated MC sim¬ 
ulations. The statistical uncertainty due to the number 
of simulated events is on the order of 0.1 %. 

The background model resulted in a scaling parame¬ 
ter of N 2 „ = 25690 for the 2v/3(3 spectrum, which is 

the best fit parameter. The uncertainty is given by the 
smallest 68 % probability interval of the marginalized 
posterior probability distribution. Using this result, the 
half-life derived according to Eq. 2 is 

T?/ 2 = (1.926 IS) • 10 21 yr ■ (3) 


6.2 Systematic Uncertainties 

The systematic uncertainties affecting the results for 
Ty, 2 were grouped into the three categories (i) detector 
parameters and fit model, (ii) MC simulation, and (in) 
data acquisition and selection. The contributions to the 
total systematic uncertainty on T 2 )/, are summarized in 
Table 2. 

(i) detector parameters and fit model 

— The systematic uncertainty on the active 76 Ge expo¬ 
sure (£av, 76 ) wa s determined using a MC approach. 
£ av. 76 is defined as 


N d ct 

£aV ,76 = Mfiif AV, if 76, i ■ (4) 

i—1 

For evaluating its uncertainty, the parameters of the 
individual detectors were randomly sampled from 
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Table 2 Contributions to the systematic uncertainty on taken into account in this work. The total systematic uncertainty 
is obtained by combining the individual contributions in quadrature. 


Item 

Uncertainty 

on Tly2 


[%] 

Active 76 Ge exposure 

±4 


Background model components 

+ 1.4 
- 1.2 


Binning 

±0.5 


Shape of the 2t'/3/3 spectrum 

< 0.1 


Subtotal fit model 


±4.3 

Precision of the Monte Carlo geometry model 

±1 


Accuracy of the Monte Carlo tracking 

±2 


Subtotal Monte Carlo simulation 


±2.2 

Data acquisition and handling 


< 0.1 

Total 


±4.8 


Gaussian distributions with mean values and stan¬ 
dard deviations according to the corresponding val¬ 
ues listed in Table 1. The correlated terms for f^v 
were also taken into account. The uncertainty on 
the live time t is 0.3%, whereas the total detec¬ 
tor masses are known with good accuracy (uncer¬ 
tainty smaller than 0.1%). The calculation yields 
£av, 76 = (13.45 ± 0.54) kg-yr. The uncertainty of 
4% is driven by the uncertainties on /a v and fjg, 
which mainly affect the number of 76 Ge nuclei in 
the active volume of the detectors, with a relatively 
smaller impact on the detection efficiency for the 
background sources. 

— The reference background model used for determin¬ 
ing accounts only for the dominant source loca¬ 
tions in the setup. The systematic uncertainty due 
to the choice of the background model components 
was evaluated by repeating the global fit with al¬ 
ternative models, which account for different source 
locations for all the background sources considered 
in the reference model. The model that accounts for 
228 Th and 228 Ac contributions also in the radon- 
shroud instead of only in the holders results in a 
1.4% longer T 2 ^. The same increase occurs if 40 K 
in the radon-shroud is added to the model compo¬ 
nents. The model including the contribution from 
214 Bi in the radon-shroud in addition to the p + sur¬ 
face and holders yields a 0.7% longer T+. In all the 
cases mentioned above, the contribution from back¬ 
ground in the 2v/3f3 spectrum region increases, since 
the peak-to-Compton ratio of the y-rays decreases 
for farther source locations leading to longer T 2 A 
estimates. Excluding contributions from very close 
source locations, like 214 Bi on the p + surface and 
60 Co on the germanium, results in a smaller increase 
of the best T 2 7 2 estimate. In this case, the contri¬ 
butions from these components are compensated by 


214 Bi and 60 Co decays in the holders, respectively. 
Consequently, the source locations are moved fur¬ 
ther out with respect to the reference model. Con¬ 
sistently, the models that include additional con¬ 
tributions from close source locations yield a de¬ 
crease in the Ty 2 value, e.g. including 214 Bi in LAr 
close to the p + surface (-1.0%) or 42 K on the ?t + 
(-1.2%) and p + (-0.6%) surfaces. Comparing alter¬ 
native background models to the reference one, the 
deviations in the Ty 2 result range between -1.2% 
and +1.4%. 

— For the standard fit, a bin width of 30 keV was used 
for the data and MC energy spectra. In order to take 
into account the systematic uncertainty related to 
binning effects, the fit was repeated twice using bin 
widths of 10 and 50 keV. The bin width of 10 keV 
was chosen in order to minimize as much as possible 
the bin size taking into account the energy resolu¬ 
tion of «4.5 keV of the coaxial detectors and the 
necessity to have enough statistics in all bins. Above 
50 keV, peak structures are washed out, leading to 
a deterioration of the fit. The deviations in the T 2 )^ 
result range between -0.5 % and +0.5 % with respect 
to that using the standard bin width. 

— The primary spectrum of the two electrons emitted 
in the 2ry3/3 decay of 76 Ge, which was then fed into 
the MC simulation, was sampled according to the 
distribution given in Ref. [33] implemented in De- 
CAyO [34]. The systematic uncertainty due to the as¬ 
sumed 2i//3/3 spectral shape was evaluated by com¬ 
paring the spectrum generated by DecayO to the 
one given in Ref. [35]. Considering the analysis win¬ 
dow used for background modeling, the maximum 
deviation is 0.2 % and the total deviation of the in¬ 
tegral in the analysis window is 0.1 %. When the fit 
with the background model is repeated using the 
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spectrum of Ref. [35] , the difference from the refer¬ 
ence If" result is less than 0.1 %. 

— A possible effect of a transition layer, where it is as¬ 
sumed that the n + dead layer on the detector sur¬ 
faces is partially active, has been investigated [36, 
37]. The dead layer thickness for individual detec¬ 
tors assumed in MC simulations were given accord¬ 
ing to the listed values in Ref. [25]. The transition 
layer is modeled using two different assumptions: a 
linearly and an exponentially increasing charge col¬ 
lection efficiency in the dead layer. The systematic 
uncertainty on T- 2 ± due to the 2i//3/3 spectrum sim¬ 
ulated with the transition layer is found to be neg¬ 
ligible. 

(ii) MC simulation 

The uncertainty related to the MC simulation arises 
from the precision of the experimental geometry model 
implemented in MaGe (1 %) and from the accuracy of 
particle tracking (2%) performed by Geant4 [38,39]. 
The total MC simulation uncertainty was estimated to 
be 2.2 % by summing in quadrature the aforementioned 
contributions. 

(in) Data acquisition and selection 
The trigger and reconstruction efficiencies for physical 
events are practically 100% above 100 keV in Gerda. 
The performance of the quality cuts applied in Phase I 
data has been investigated through a visual analysis. 
The total uncertainty related to data acquisition and 
selection was estimated to be less than 0.1 %. 

Summing in quadrature the uncertainties of the three 
groups gives a total systematic uncertainty of ±4.8 %. 


6.3 Results and Discussion 


Fig. 2 shows the experimental data together with the 
best fit model for the golden data set. The different 
components of the minimum background model are also 
reported. The model is able to reproduce the experi¬ 
mental data well, as shown in the lower panel of the 
figure by the residuals. 

The best estimate of the Ty. 2 of the 2is/3/3 decay of 
76 Ge is: 


^1/2 (1.926 io.022 stat -0.092 systj 

= (1.926 ±0.095) • 10 21 yr, 


10 21 yr 


(5) 


with the latter combining in quadrature the statistical 
(fit) and systematic uncertainties. The total uncertainty 
of 4.9 % is dominated by the systematic uncertainties. 
The largest contribution to the systematic uncertainties 
comes from the uncertainty on the active 76 Ge exposure 


(4%), which can only be reduced by performing new 
and more precise measurements of the active masses 
of the coaxial detectors. Other significant contributions 
are related to the Monte Carlo simulations (2.2 %) and 
to the background model assumptions (+i' 2 %)- The lat¬ 
ter have been significantly reduced in this analysis com¬ 
pared to the analysis of the first 5 kg-yr of Phase I data 
reported in Ref. [41], where the systematic uncertainty 
due to the background model was The new re¬ 

sult is in good agreement with that mentioned above. 
Adding further identified components to the reference 
background model results in a slight increase of the best 
Ty, 2 estimate. 

The background level achieved in Gerda Phase I 
is about one order of magnitude lower with respect 
to predecessor ,6 Ge experiments, and has allowed the 
measurement of Tj 2 7 2 with an unprecedented signal-to- 
background ratio of 3:1 in the 570-2039 keV interval. 
The ratio amounts to 4:1 for the smaller interval of 
600-1800 keV. 


7 Limits on Majoron-emitting double /3-decays 
of 76 Ge 

7.1 Analysis 

The search for Ovfifdx was performed using the golden 
and BEGe data sets, amounting to a total exposure 
of 20.3 kg-yr. The analysis employed the background 
model described in section 5. The information from the 
two data sets was combined in one fit, while keeping 
their energy spectra distinct. A separate fit was per¬ 
formed for each spectral index, containing the back¬ 
ground contributions, the contributions from 2i/f3j3 de¬ 
cay, and also the Majoron component under study. A 
single parameter, Ty.f ■ is considered common for the 
two data sets. It is defined as the half-life of the respec¬ 
tive Majoron accompanied mode. 

In order to improve the detection efficiency for the 
Majoron processes with low n (n = 1,2), a slightly 
different event selection was used with respect to the 
T 2 ± analysis. If an event occurs with energy deposition 
in two detectors and the energy deposit in the detec¬ 
tor where the decay took place is below the threshold 
for the anti-coincidence cut, the event contributes to 
the energy spectrum of the other detector. Therefore, 
when determining the total energy spectrum resulting 
from decays in one of the detectors, the energy spectra 
from all detectors in the array have to be taken into 
account. Such a selection has no impact on the detec¬ 
tion efficiency for the Majoron process with n = 3 and 
7 and 2i//3/3 decay. The content of the i-th bin in the 
combined energy spectrum of all N ( j et detectors in the 
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1000 1500 2000 

energy (keV) 


Fig. 2 Upper panel: experimental data (markers) and the best fit model (black histogram) for the golden data set. The 
contribution from 2u/3/3 (green) and from the single background components are also shown. Lower panel: ratio between 
experimental data and the prediction of the best fit model. The green, yellow and red regions are the smallest intervals 
containing 68%, 95% and 99% probability for the ratio assuming the best fit parameters, respectively [40]. 


array, for decays taking place in the active and dead 
part of detector a , becomes: 




(In 2) N a 


m e 


■T 


'OVX 


M a /76,o 


1/2 


N det 

fAV,a tj £ AV,j@ 

j =1 


a,0vx 

AV,i,j 


N det 

+ (1 — fAV,a) ^2, tj £ DL,j& 


3 =1 


a, Ovx 
DL,i,j 


( 6 ) 


with ' ) °Dhfj ) giving the content of the i-th bin 

of the normalized energy distribution recorded with de¬ 
tector j for Qvppx taking place in the active (dead) 
volume of detector a. Summing up the simulations of 
decays in all Ndet detectors results in the final model 


spectrum: 


N ds 


=£ a- 


Oz'X 


a—1 


(7) 


For all four Majoron modes (n = 1,2,3, 7) only lower 
limits on the half-life can be given. They were obtained 
from the 90 % quantiles of the marginalized posterior 
distributions. These lower limits for n °f taking 

into account the systematic uncertainties, are in units 
of 10 23 yr: >4.4, >1.9, >0.9, and >0.4 for n = 1, 2, 3, 
and 7, respectively. The respective half-life of the 2i//3/3 
process derived from this analysis amounts to in units of 
10 21 yr: 1.96±0.03 stat , 1.97±0.03 stat , 1.98±0.03 stat , and 
1.99±0.03 s t a t- Within the uncertainties coming from 
the different background models and the different data 
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sets of the two analyses, the derived T(J. 2 values are in 
agreement (<lcx) with that discussed in section 6.3. 

7.2 Systematic Uncertainties 

The systematic uncertainties were divided into the three 
categories (i) detector parameters and fit model , (ii) MC 
simulation , and (in) data acquisition and selection. 

(i) detector parameters and fit model 
Uncertainties from the fitting procedure were folded 
into the posterior distribution of with a MC ap¬ 
proach. Each source of uncertainty is described by a 
probability distribution. The fitting procedure was re¬ 
peated 1000 times, each time drawing a random number 
for each source of uncertainty according to its proba¬ 
bility distribution: 

— Material screening measurement results were used 
to constrain the minimum number of events expected 
from close and medium distant sources of the 214 Bi 
and 228 Th decays. Gaussian distributions describing 
these lower limits used in the fit were derived from 
the mean and standard deviations of the screening 
measurements. For details refer to Ref. [25]. 

— As for the T(Y 2 analysis, the standard fit uses a bin 
width of 30 keV for the data and MC energy spectra. 
In order to determine the systematic uncertainty re¬ 
lated to binning effects the bin width was sampled 
uniformly from 10 keV to 50 keV. 

— Uncertainties on the active volume fractions enter 
the model in several ways. On the one hand, the 
MC energy spectra for all internal sources, that is 
for 2v(3/3, 0i//3/?x, 60 Co, and 68 Ga decays, are af¬ 
fected, as the fraction of decays taking place in the 
active and dead part of the detectors changes with 
changing f.\v- On the other hand, the uncertainty 
on the active volume fraction also plays a role for 
the shape of the energy spectrum due to 42 K decays 
on the n + surface. Larger /av means thinner n + 
dead layer and thus the possibility of an increased 
contribution from the electrons to the spectrum. For 
smaller f av and thicker n + dead layer, their con¬ 
tributions are expected to be reduced. The active 
volume fraction for each detector was sampled from 
a Gaussian distribution with mean and standard de¬ 
viation according to Table 1. For the coaxial detec¬ 
tors, the partial correlations of the uncertainty were 
taken into account. The simulated spectra of the in¬ 
ternal sources as well as of the 42 K decays on the 
n + surface are composed according to the sampled 
active volume fractions. 

— The uncertainty on the fraction of enrichment in 
76 Ge of the germanium that constitutes the detec¬ 


tors plays a role when converting the number of 
events attributed to 0^/3/3x into The prob¬ 

ability distribution of fi§ for each detector is given 
by a Gaussian function with mean values and stan¬ 
dard deviations as listed in Table 1. 

— The data does not allow the resolution of the ambi¬ 
guity regarding the exact positions of the near and 
medium distant sources. The 214 Bi decays serves as 
a representative in order to estimate the impact of 
this uncertainty. Their near position is represented 
by decays in the holders, in the mini-shroud or on 
the n + surface of the detectors, each having a prob¬ 
ability of 1/3 in the sampling process. The medium 
distant position is represented by decays in the radon- 
shroud or in the LAr, having a probability 1/2 in 
contrast. 

— Extensive studies of the characteristics of the BEGe 
diodes suggest the presence of a transition layer be¬ 
tween the region where the detector is fully efficient 
and the external dead region [36,37]. An uncertainty 
as high as ±0.5% on the lower limits of is es¬ 
timated for this effect in the case of the BEGe de¬ 
tectors. This uncertainty was folded into the total 
marginalized posterior distribution a posteriori. The 
corresponding uncertainty for the coaxial detectors 
is estimated to be negligible. 

The marginalized posterior distributions for Tffff 
derived from each of the 1000 individual fits were sum¬ 
med up. The resulting total marginalized posterior dis¬ 
tribution accounts for the statistical as well as for the 
listed systematic uncertainties related to the fit model. 

As for the T// 2 analysis, the uncertainties on the ac¬ 
tive volume fractions and on the enrichment fractions 
are major contributions to the total uncertainty on the 
limits for Ty(f. However, the largest source of uncer¬ 
tainty is the composition of the fit model and the indi¬ 
vidual background contributions. In the case of n = 1, 
a fit with a bin width of 50 keV weakens the limit by 
~ 16 % compared to the standard fit, while the result 
for T 2 / 2 is not affected at all. The stability of the T 2 U 
results shows the validity of the fit. The use of the al¬ 
ternative close and medium distant source positions for 
214 Bi decays leads to maximal variations of /® 2 3 6 % °f 
the limit on T^ 2 X . 

(ii) MC simulation 

As in the case of the Xj 2 / 2 measurement, a total MC 
simulation uncertainty of 2.2% has to be taken into 
account for effects related to the geometry implemen¬ 
tation and particle tracking. It is folded into the total 
marginalized posterior distributions. No effect on the 
lower limits is observed for any of the spectral modes. 

(in) Data acquisition and selection 
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Table 3 Experimental results for the limits on of 76 Ge for the Majoron models given in Refs. [7,42,43,44], The first 

section considers lepton number violating models (I) allowing Ou/3/3 decay, while in the second section lepton number conserving 
models (II) are listed, where Oi'/3/3 decay is not allowed. The first column gives the model name, the second the spectral index, 
n, the third the information on whether one Majoron, Xi or two Majorons, XX-* is emitted, the fourth if the Majoron is a 
Goldstone boson, the fifth provides its lepton number, L, the sixth the experimental limit on T 9 J 2 X of 76 Ge obtained in this 
analysis. The nuclear matrix elements, A4 0ux , the phase space factor, G°^ x , and the resulting effective coupling constants, 
(g), are given in the seventh, eighth and ninth columns, respectively. The limits on of 70 Ge for the Majoron models and 

{g) correspond to the 90% quantiles of the marginalized posterior probability distribution. For the case of n = 1, the nuclear 
matrix element, _A/| 0l ' x , from Refs. [45,46,47,48,49,50,51] and the phase space factor, G 0l ' x , from Ref. [52] are used for the 
calculation of ( g ). The given range covers the variations of M 0vx in these works. For n = 3 and 7, {g) is determined using 
the matrix elements and phase space factors from Ref. [42]. The results for OvPPx {n = 3, 7) account for the uncertainty on 
A4 0ux . For n = 2, only the experimental upper limit is given. 


Model 

n 

Mode 

Goldstone 

boson 

L 

rjiOb'X 

i l/2 

[10 23 yr] 

M 0ux 

qOvx 

[yr- 1 ] 

(9) 

IB 

IC 

ID 

IE 

IF 

i 

1 

3 

3 

2 

X 

X 

XX 

XX 

X 

no 

yes 

no 

yes 

bulk field 

0 

0 

0 

0 

0 

> 4.2 

> 4.2 

> 0.8 

> 0.8 

> 1.8 

(2.30 - 5.82) 
(2.30 - 5.82) 
10“ 3±1 
10“ 3±1 

5.86- 10" 17 
5.86- 10" 17 
6.32 • 10" 19 
6.32 • 10" 19 

< (3.4 - 8.7) • 10“ 5 

< (3.4-8.7) • 10 -5 

< 2.1± 4 '4 

< 2.ilti 

IIB 

1 

X 

no 

-2 

> 4.2 

(2.30 - 5.82) 

5.86 • 10" 17 

< (3.4-8.7) • 10 -5 

IIC 

3 

X 

yes 

-2 

> 0.8 

0.16 

2.07- 10" 19 

< 4.7- 10" 2 

IID 

3 

XX 

no 

-1 

> 0.8 

10“ 3±1 

6.32 • 10" 19 

< 2.1 

HE 

7 

XX 

yes 

-1 

> 0.3 

10“ 3±1 

1.21 • 10" 18 

< 2.2± 4 '4 

IIF 

3 

X 

gauge boson 

-2 

> 0.8 

0.16 

2.07- 10" 19 

< 4.7- 10" 2 


The uncertainty from data acquisition and selection is 
estimated to be below 0.1 % and does not alter the de¬ 
rived limits on T®J X . 

7.3 Results and Discussion 

Fig. 3 shows the global model for the case of spectral in¬ 
dex n = 1 together with the energy spectra for both the 
coaxial and the BEGe data sets. The contributions from 
the background contaminations, from the 2vf3j3 decay 
only, and the combined spectra from the background 
contaminations and 2v/3/3 decay are drawn separately. 
The 35868 events in the data spectrum of the golden 
data set were matched with 35834 events in the best- 
fit model for n = 1. Of those events, in the best fit, 
54.5 are attributed to Ovfipx- For the BEGe data set, 
the best-fit model contains 5081.4 counts for the 5035 
measured events. In this fit, 7.8 events are attributed 
to OvpPx decay. The limit of T^ x at 90% C.I. derived 
from the fit is also drawn (green histogram). The up¬ 
per limits at 90% C.I. for the remaining three modes 
are reported for illustrative purpose (blue histogram for 
n = 2, orange for n = 3 and red for n = 7). The maxi¬ 
mum of the corresponding distributions shifts to higher 
energy with the diminishing of the spectral index n. 
The resulting lower limits on T^j x , determined as the 
90 % quantiles of the posterior probability distributions 
and taking into account all uncertainties related to the 


fit model, are (in units of 10 23 yr): >4.2, >1.8, >0.8 
and >0.3 for n = 1,2,3 and 7, respectively. The results 
are summarized in Table 3 for the different Majoron 
models. 

The limits on T®J X presented here are the most 
stringent limits obtained to date for ,6 Ge. The limits 
for n = 1 and n = 3 are improved by more than a fac¬ 
tor six [9], the limit for n = 7 is improved by a factor 
five [8] compared to previous measurements. The limit 
for the mode with n = 2 is reported here for the first 
time. 

From the lower limits on T®J X , upper limits on the 
effective neutrino-Majoron coupling constants (g) for 
the models with n = 1,3 and 7 can be calculated using 
the following equations: 

1/7?;* = |( 5 )| 2 • G°" X (Q 00 , Z) ■ \M 0ux \ 2 (8) 

and 

1/T?;* = |( 5 )| 4 • G Ovxx (Q 00 , Z) ■ \M 0uxx \ 2 (9) 

for single and double Majoron emission, respectively. 
The matrix element for the models with n = 1 (IB, 
IC and IIB) are taken from Refs. [45,46,47,48,49,50, 
51], whereas the phase space factor is that of Ref. [52], 
The matrix elements for the models with n = 3 (ID, 
IE, IIC, IID, IIF) and with n = 7 (HE) as well as 
the corresponding phase space factors are taken from 
Ref. [42] . The results for the upper limits on ( g) are also 
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Fig. 3 Best-fit model and data energy spectrum for the coaxial and the BEGe data sets for the case of spectral index 
n — 1. The contributions from 2u/3/3 decay and the background contributions are shown separately. The best-fit model does 
not contain the contributions from 0v/3/3x- The smallest interval of 68 % probability for the model expectation is indicated in 
grey. Also shown is the upper limit for 0u/3/3x decay with n = 1 as determined from the 90 % quantile of the marginalized 
posterior probability for F° r illustrative purpose, also the upper limits at 90% C.l. of the other three spectral indices 

n = 2, 3, 7 are reported. 


shown in Table 3. The coupling constants allow a com¬ 
parison with other isotopes. The best limits on OvftPx 
decay of isotopes other than 76 Ge have been obtained 
for 100 Mo [10] and 136 Xe [14]. When comparing with 
the case of 100 Mo, it becomes obvious that the limits 
on TfT* determined in the present analysis are about 
one order of magnitude more stringent, for the case of 
n = 7 even two orders of magnitude. However, due to 
the differences in the matrix elements and the phase 
space factors, the resulting limits on ( g ) from lon Mo 
and ,6 Ge are comparable. The limits for ( g ) derived 
from 136 Xe are a factor of two to five more stringent 


due to the higher limits that had been measured for 

1 i/2 ■ 

8 Conclusions 

Phase I of the Gerda experiment, located at the INFN 
Laboratori Nazionali del Gran Sasso (LNGS) in Italy, 
has been executed between November 2011 and May 
2013. Utilizing the collected exposure of Phase I, an im¬ 
proved result of the half-life of the 2v/3P process in ,6 Ge 
was obtained and new limits for the half-lives of the 
Majoron-emitting double beta decays were produced. 
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The half-life for the 2i//3/3 process is determined to 
be: 

Ty 2 = (1.926 ± 0.095) • 10 21 yr . (10) 

Thanks to the extremely low background level in the 
Gerda experiment, with a signal-to-background ratio 
of 3:1 in the 570-2039 keV interval and a refined back¬ 
ground model, the measurement has an unprecedented 
precision (<5%) with respect to previous experiments 
using ,6 Ge. The new result is in good agreement with 
the one derived from a smaller data set with 5 kg-yr ex¬ 
posure [41]. The inclusion of more components into the 
reference background model results in a slight increase 
of the best estimate for Ty 2 ■ 

Majoron emission processes were searched for in the 
energy spectra using an exposure of 20.3 kg-yr. The 
analysis was performed for all four possibilities of the 
spectral index n {n = 1, 2, 3, and 7). No indication 
for a contribution of Ov/lfix was found in any of the 
cases. Lower limits on the half-lives, Ty 2 , were deter¬ 
mined from the quantiles of 90 % probability of the 
marginalized posterior probability distributions. The 
results constitute the most stringent limits on of 
76 Ge obtained to date. For the standard mode (n = 1), 
the lower limit is determined to be: 

T°J* > 4.2 • 10 23 yr. (11) 

From the lower limit on T-Jr*, an upper limit on the 
effective neutrino-Majoron coupling constant, ( g ), can 
be inferred: 

(g) < (3.4- 8.7) • 1(T 5 . (12) 
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